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Abstract: Near critical catalyst-reactant branching processes with controlled immigra- 
tion are studied. The reactant population evolves according to a branching process whose 
branching rate is proportional to the total mass of the catalyst. The bulk catalyst evolu- 
tion is that of a classical continuous time branching process; in addition there is a specific 
r^ ' form of immigration. Immigration takes place exactly when the catalyst population falls 

' *~i , below a certain threshold, in which case the population is instantaneously replenished to 

the threshold. Such models are motivated by problems in chemical kinetics where one 
wants to keep the level of a catalyst above a certain threshold in order to maintain a 
desired level of reaction activity. A diffusion limit theorem for the scaled processes is 
presented, in which the catalyst limit is described through a reflected diffusion, while the 
^ \ reactant limit is a diffusion with coefficients that are functions of both the reactant and 

^\ • the catalyst. Stochastic averaging principles under fast catalyst dynamics are established. 

[■*-■ ' In the case where the catalyst evolves "much faster" than the reactant, a scaling limit, 

OO , in which the reactant is described through a one dimensional SDE with coefficients de- 

>0 ' pending on the invariant distribution of the reflected diffusion, is obtained. Proofs rely 

rf) ' on constrained martingale problem characterizations, Lyapunov function constructions, 

(^ , moment estimates that are uniform in time and the scaling parameter, and occupation 

^Nj ' measure techniques. 
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1. Introduction 



This work is concerned with catalytic branching processes that model the dynamics of catalyst- 
reactant populations in which the activity level of the reactant depends on the amount of 
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catalyst present. Branching processes in catalytic environments have been studied extensively 
and are motivated, for instance, by biochemical reaction networks (see [4, 6, 13] and references 
therein). The basic model consists of a continuous time countable state branching process de- 
scribing the evolution of the catalyst population and a second branching process for which the 
branching rate is proportional to the total mass of the catalyst population modeling the evolu- 
tion of reactant particles. Such processes were introduced in [4] in the setting of super-Brownian 
motions ([13]). For classical catalyst-reactant branching processes, the catalyst population dies 
out with positive probability and subsequent to the catalyst extinction, the reactant popula- 
tion stays unchanged and therefore the population dynamics are modeled until the time the 
catalyst becomes extinct. In this work, we consider a setting where the catalyst population 
is maintained above a positive threshold through a specific form of controlled immigration. 
Branching process models with immigration have also been well studied in literature (see 
[1, 13] and references therein). However, typical mechanisms that have been considered corre- 
spond to adding an independent Poisson component (see e.g. [10]). Here, instead, we consider a 
model where immigration is allowed only when the population drops below a certain threshold. 
Roughly speaking, we consider a sequence {X(")}„gN of continuous time branching processes, 
where X'"' starts with n particles. When the population drops below n, it is instantaneously 
restored to the level n. 

There are many settings where controlled immigration models of the above form arise nat- 
urally. One class of examples arise from predator-prey models in ecology, where one may be 
concerned with the restoration of populations that are close to extinction by reintroducing 
species when they fall below a certain threshold. In our work, the motivation for the study of 
such controlled immigration models comes from problems in chemical kinetics where one wants 
to keep the level of a catalyst above a certain threshold in order to maintain a desired level of 
reaction activity. 

Our main goal is to establish diffusion approximations for such catalyst-reactant systems 
under suitable scaling. We consider two different scaling regimes; in the first setting the catalyst 
and reactant evolve on "comparable timescales" , while in the second setting the catalyst evolves 
"much faster" than the reactant. In the former setting, the limit model is described through a 
coupled system of reflected stochastic differential equations with reflection in the space [1, oo) x 
M. The precise result (Theorem 2.1) is stated in Section 2. For the latter setting, where the 
catalyst evolution is much faster, we establish a stochastic averaging limit theorem. A key 
ingredient here is an ergodicity result, which says that under a suitable "criticality from below" 
assumption on the catalyst dynamics, the limiting catalyst reflected diffusion admits a unique 
stationary distribution, which takes an explicit form (Proposition 3.1). Characterization of the 
invariant distribution is based on a variant of Echeverria's criterion for constrained Markov 
processes [12]. Next, by constructing suitable uniform Lyapunov functions, we show that the 
stationary distribution of the scaled catalyst branching process converges to that of the catalyst 
diffusion (Theorem 3.1). These results are then used to establish a stochastic averaging principle 
that governs the dynamics of the reactant population in the fast catalyst limit. The limit 
evolution of the reactant population is given through an autonomous one dimensional SDE 
with coefficients that depend on the stationary distribution of a reflected diffusion in [l,oo). 
Such model reductions are important in that they not only help in better understanding the 
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dynamics of the system but also help in reducing computational costs in simulations. We refer 
the reader to [9] and references therein for similar results in the setting of chemical reaction 
networks. Proofs proceed by developing suitable moment estimates that are uniform in time 
and the scaling parameter and by using characterization results for probability laws of reflected 
diffusions through certain constrained martingale problems [11]. 

We consider two different formulations of models with multiple time scales. In Theorem 
4.1 we consider the setting where both catalyst and reactant processes are described through 
(reflected) diffusions and the time scale parameter appears in the coefficients of the catalyst evo- 
lution equation. An important step here is to argue that the generator of the two dimensional 
catalyst-reactant branching process is suitably close to the generator of the one dimensional 
averaged diffusion, for large values of the scaling parameter. Bounds on the exponential mo- 
ments of the catalyst process, obtained in Lemma 8.1, play a key role in this argument. The 
second formulation is considered in Theorem 4.2. Here, both catalyst and reactant populations 
evolve according to near critical countable state branching processes and the branching rate 
in the catalyst dynamics is of higher order than that for the reactant process. In this setting 
one encounters the additional difficulty of showing that the steady state distributions of the 
scaled catalyst branching process, for large values of the scaling parameter, are suitably close 
to the stationary distribution of the limiting catalyst reflected diffusion. The approach taken 
here is based on characterizing the limit points of a certain sequence of random measures on 
the path space of the catalyst process and the associated reflection process, as time and the 
scaling parameter together approach infinity. 

The paper is organized as follows. We begin in Section 2 by presenting the basic limit theorem 
in the setting of "comparable time scales" . Section 3 studies the time asymptotic behavior of 
the catalyst process under a suitable criticality from below assumption. Section 4 presents our 
main results for the multiple time scale setting. Section 5 collects some auxiliary estimates 
that are needed in our proofs. Section 6 proves Theorem 2.1 and Section 7 is devoted to the 
proofs of Proposition 3.1 and Theorem 3.1. Finally, in Section 8 we present proofs of stochastic 
averaging principles stated in Section 4. 



1.1. Notation 



The following notation will be used throughout this work. Denote by N the natural numbers, let 
No := Nu{0}, denote the set of integers by Z, and let ]R+ := [0, oo) be the set of nonnegative real 
numbers. The state spaces of the scaled catalyst, reactant, and auxiliary processes, X^"'',Y^"'', 
and Z("), respectively, introduced below in (2.1), are S^^ := {^|/ G No} n [l,oo), E^ := 

{^\l G No}, and E^ := {^\l G Z}. Let W^ := sj^ x sj"^ x S^"^ and W := [1, oo) x M+ x M. 
Let C'^(W) denote the space of A;-times continuously differentiable, real valued functions on 
W, and denote by C^(W) the space of C (W) functions with compact support. Here, by a 
(/c-times) differentiable function / on a set D C M" we mean a function that can be extended 
to a (/s-times) differentiable function / on an open domain U D D such that / restricted to 
D equals /. Given a metric space S, the space of probability measures on S will be denoted 
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by V{S), the Borel cr-field on 5 by 13{S), and the space of real valued, bounded, measurable 
functions on S by BM(S'). Let 

-D(M4. : 5) := {/ : M+ — )• 5|/ is right continuous and has left limits} 

and Di{R+ : M) := {/ G L>(M+ : M)|/(0) > 1}, where these D-spaces are endowed with the 
usual Skorohod topology. Let C(M+ : M+) be the space of continuous functions from M+ to M+ 
endowed with the local uniform topology. We say a sequence {^njnsN of random variables with 
values in some Polish space £ is tight if the corresponding probability laws are a tight sequence 
in V{£). For a function ^ : M_|_ — )• M", let the jump at time t be defined as A^j := ^t — £,t-, 
t > 0, and A^o := 0. For a function / : M+ -^ R and t > 0, let |/|*,j := sups<j |/(s)|. For two 
semimartingales ^ and C, the quadratic covariation (or bracket process) and predictable (or 
conditional) quadratic covariation are denoted by {['?,C]t}i6K+ ^nd {(^,C)t}t6K+ respectively; 
their definition will be recalled in Section 5. 



2. Diffusion Limit under Comparable Timescales 



Consider a sequence of pairs of continuous time, countable state Markov branching processes 
(X^"\y("^), where X^^' and y(") represent the number of catalyst and reactant particles, 

(n) 

respectively. The dynamics are described as follows. Each of the X^ particles alive at time t 
has an exponentially distributed life time with parameter A^" (mean life time 1/A^" ). When it 
dies, each such particle gives rise to a number of offspring, according to the offspring distribution 
//-["(•). Additionally, if the catalyst population drops below n, it is instantaneously replenished 
back to the level n {controlled immigration). The branching rate of the reactant process Y^"^' 
is of the order of the current total mass of the catalyst population, i.e. X^'^> /n, and we denote 
the offspring distribution of Y^^' by (jl^P {■). A precise definition of the pair (X^^' ,Y^^') will 
be given below. We are interested in the study of asymptotic behavior of {X^"-> ,Y^"->), under 
suitable scaling, as n — )• oo. 

To facilitate some weak convergence arguments, we will consider an auxiliary sequence of 
processes Z'"^ that "shadow" X^"^' in the following manner. The process Z^"'> will be a Z 
valued pure jump process whose jump instances and sizes are the same as that of X '"^ away 
from the boundary {n}, whereas when X^"' is at the boundary, Z^"'' has a negative jump of 
size 1 whenever there is immigration of a catalyst particle into the system. This description 
is made precise through the infinitesimal generator given in (2.2). The process Z^"^' will not 
appear in the statements of the results, nevertheless it plays an important role in our proofs. 

We now give a precise description of the various processes and the scaling that is considered. 
Roughly speaking, time is accelerated by a factor of n and mass is scaled down by a factor of 
n. Define RCLL processes 



'{'") yin) r,{ri) - 
'nt nt nt 

n ' n ^ n 



^{n) _^ (x^n)^^{n)^^(n)\ _^ ( A^^ £n^^ £ni_ ) , t € M+, (2.1) 
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and let Wq = (xq , t/q , Zq ) G W^"\ where {nx^ ^ny^ ) is the initial number of catalyst 
and reactant particles and Zq = Xq . Then {Wj }tgK_^ is characterized as the W'-"-* valued 

Markov process with sample paths in Z)(]R_|. : W'"'), starting at Wq "^ = (xg , ^q , Zq ), and 
having infinitesimal generator A}"^' given as 



fc=0 

oo u ^ 

+ x'^^\^xyY,[t>{x,y + —-,z)-ct>{^)\lif\k), (2.2) 



fc=0 



where w = {x,y,z) G W^"^ and <p G BM(W). From the definition of the generator we see 
that, for each k > 0, given Wj" = {x,y,z) G W^"^ the process jumps to {x,y + ^^^,z) 
with rate X2 n'^xfi2 {^)i ^^'^ to {x + -^,y,z -\ — ^) with rate A^" n^/i]^" (/c), except when 
k = and x = 1, in which case the latter jump is to {x,y,z + — ^) with rate Xl n^/i^ (0). 
This property of the generator at x = 1 accounts for the instantaneous replenishment of the 
(unsealed) catalyst population to level n, whenever the catalyst drops below n. 

For i = 1,2, let 

00 00 00 

mt^ ■■= E ^/^^^(^)' "S"^ = E(^- - i)vr^(A^), 4^ = E(^ - -!"^)V?^(A.). 

fc=0 fc=0 fc=0 

We make the following basic assumption on the parameters of the branching rates and offspring 
distributions as well as on the initial configurations of the catalyst and reactant populations. 

Condition 2.1. (i) Fori =1,2 and for n G N, a^\x^ G (0, 00) andrn^' = 1+^, cf G 
(— n, 00). 
(ii) For i = 1,2, as n ^ 00, c^ — )• Cj G M, a^ — )■ a^ G (0, 00), and A^ — )■ Aj G (0, 00). 
(in) For i = 1,2 and for every e G (0, 00), lim„^oo Y^uye^i^ ~ "^i I'k^i V') = 0- 
(iv) As n^ 00, {x}^\yo) -^ ixo,yo) G [l,oo) x M+. 

Condition 2.1 and the form of the generator in (2.2) ensure that the scaled catalyst and 
reactant processes transition on comparable time scales, namely O(n^). In order to state the 
limit theorem for {X^^' jY^"^'), we need some notation and definitions associated with the one 
dimensional Skorohod map with reflection at 1. Let F : Z)i(]R+ : M) — t- 15(1^+ : [l,oo)) be 
defined as 

F(V')(t) := {^p{t) + 1) - inf {V'(s) A 1}, for V G D{R+ : E). (2.3) 

0<s<t 

The function F, known as Skorohod map, can be characterized as follows: If ip, (p, Vj* G L'(M+ : 
M) are such that (i) ^(0) > 1, (ii) (j) = ip + Vj* , (iii) </> > 1, (iv) rf is non-decreasing, 
/fOoo) -'-{0(5)^1} '^^*(^) ~ *-*' ^^"^ ^*(0) = 0' then (j) = F(V') and t]* = (p — ip. The process r]* 
can be regarded as the reflection term that is applied to the original trajectory ip to produce a 
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trajectory (f) that is constrained to [1, oo). We will make use of the following Lipschitz property 
of the Skorohod map: For ip,ijj £ D{M.^ : M), 

sup|r(V')(s) -r(^)(s)| < 2sup|^(s) -^(s)|. (2.4) 

s<t s<t 



The diffusion limit of (X^^^y^")) will be the process {X,Y), starting at {xo,yo), which is 
given through a system of stochastic integral equations as in the following proposition. 

Proposition 2.1. Let {Ct,^,P,{^t}) be a filtered probability space on which are given in- 
dependent standard {^t} Brownian motions B^ and B^ . Let Xq,Yq be square integrable J-q 
measurable random variables with values in [l,oo) and M+ respectively. Then the following 
system of stochastic integral equations has a unique strong solution: 

Xt = r(xo+ f ciXiXsds + / y^aiXiXsdBA (t), (2.5) 

Yt = Yo+ [ C2\2XsYsds + [ ^a2X2XsYsdBj, (2.6) 

Jo Jo 

r]t = Xt-Xo- [ ciXiXsds - I ^aiXiXsdBf, (2.7) 

Jo Jo 

where T is the Skorohod map defined in (2.3). 



In the above proposition, by a strong solution of (2.5) - (2.7) we mean an ^-adapted 
continuous process {X,Y,rj) with values in [l,oo) x M_|_ x M4. that satisfies (2.5) ~ (2.7). The 
following is the main result of this section. 

Theorem 2.1. Suppose Condition 2.1 holds. Then the process {X^"^' ,Y^'^') converges weakly 
in D(M+ : [l,oo) x M+) to the process {X,Y) given in Proposition 2.1 with {Xq,Yq) = (a;o,yo)- 

Proposition 2.1 follows by standard arguments, so its proof is relegated to the appendix. 
Theorem 2.1 will be proved in Section 6. 



3. Asymptotic Behavior of the Catalyst Population 



Stochastic averaging results in this work rely on understanding the time asymptotic behavior 
of the catalyst process. Such behavior, of course, is also of independent interest. We begin 
with the following result on the stationary distribution of X, where X is the reflected diffusion 
from Proposition 2.1, approximating the catalyst dynamics (Theorem 2.1). The proof uses an 
extension of the Echeverria criterion for stationary distributions of diffusions to the setting 
of constrained diffusions (see Section 7.1). We will make the following additional assumption. 
Recall the constants c\ G (— n,cx)) and ci € M introduced in Condition 2.1. 

Condition 3.1. For alln £^, c\ ' <{) and ci < 0. 
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Proposition 3.1. Suppose Condition 3.1 holds. The process X defined through (2.5) has a 
unique stationary distribution, ui, which has density 

,(.).J|-pPS-). ■^-•SJ (3,, 

0, ifx<l, 



where 9 := (^f^{lexpi2^^x))dx^ \ 



The following result shows that the time asymptotic behavior of the catalyst population is 
well approximated by that of its diffusion approximation given through (2.5). We make the 
following additional assumption on the moment generating function of the offspring distribu- 
tion, which will allow us to construct certain "uniform Lyapunov functions" (see Theorem 7.2 
and the function V^"^' defined in (7.5)) that play a key role in the analysis. 

Condition 3.2. For some 6 > 0, 

oo 

supye^'Vi"^(A;)<oo. (3.2) 

Theorem 3.1. Suppose Conditions 2.1, 3.1, and 3.2 hold. Then, for each n G N, the process 
X^"^' has a unique stationary distribution Vi , and the family {z^}"" }neN is tight. As n ^ oo, 
z^l" converges weakly to vi. 

Proposition 3.1 and Theorem 3.1 will be proved in Section 7. 



4. Diffusion Limit of the Reactant under Fast Catalyst Dynamics 

As noted in Section 2, the catalyst and reactant populations whose scaled evolution is described 
through (2.2) transition on comparable time scales. In situations in which the catalyst evolves 
"much faster" than the reactant, one can hope to find a simplified model that captures the 
dynamics of the reactant population in a more economical fashion. One would expect that 
the reactant population can be approximated by a diffusion whose coefficients depend on the 
catalyst only through the catalyst's stationary distribution. Indeed, we will show that the 
(scaled) reactant population can be approximated by the solution of 

Yt = Yo+ C2\2mxYsds+ J a2X2mxYsdBs, Yq = yo, (4.1) 

Jo Jo 

where nix = J^ xui{dx) = — ^^ exp(2ci/ai). 

Such model reductions (see [9] and references therein for the setting of chemical reaction 
networks) not only help in better understanding the dynamics of the system but also help in 
reducing computational costs in simulations. In this section we will consider such stochastic 



/Near Critical Catalyst Reactant Branching Processes with Controlled Immigration 8 

averaging results in two model settings. First, in Section 4.1, we consider the simpler setting 
where the population mass evolutions are described through (reflected) stochastic integral 
equations and a scaling parameter in the coefficients of the model distinguishes the time scales 
of the two processes. In Section 4.2 we will consider a setting which captures the underlying 
physical dynamics more accurately in the sense that the mass processes are described in terms 
of continuous time branching processes, rather than diffusions. 



4^.1. Stochastic Averaging in a Diffusion Setting 



In this section we consider the setting where the catalyst and reactant populations evolve 
according to (reflected) diffusions similar to X and Y from Proposition 2.1, but where the 
evolution of the catalyst is accelerated by a factor of n (i.e. drift and diffusion coefficients 
depend on n). More precisely, we consider a system of catalyst and reactant populations that 
are given as solutions of the following system of stochastic integral equations: Fot t > 0, 

Xf ) = r Ix'^^^ + f nciXiXi''Us+ f \l naiXiXt'^ dsA (t) 

Jo Jo 

where {Xq, Iq ) = (xq, yo), ci, C2 S ffi, Oj, Aj G (0, oo), B^ and B^ are independent standard 
Brownian motions, and T is the Skorohod map described above Proposition 2.1. 

The following result says that if ci < 0, then the reactant population process Y^"'' , which is 
given through a coupled two dimensional system, can be well approximated by the one dimen- 
sional diffusion Y in (4.1), whose coefficients are given in terms of the stationary distribution 
of the catalyst process. 

Theorem 4.1. Suppose Condition 3.1 holds. The process Y^"^' converges weakly in ' 
to the process Y . 

The proof of Theorem 4.1 is given in Section 8. 



4.2. Stochastic Averaging for Scaled Branching Processes 



We now consider stochastic averaging for the setting where the catalyst and reactant popula- 
tions are described through branching processes. Consider catalyst and reactant populations 
evolving according to the branching processes introduced in Section 2, but where the catalyst 
evolution is sped up by a factor of n. That is, we consider a sequence of catalyst populations 
X^ := X!^ , t > 0, where X^"'^ are the branching processes introduced in Section 2. The 
reactant population evolves according to a branching process, Y^^\ whose branching rate. 
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as before, is of the order of the current total mass of the catalyst population, X'"-*/^- The 
infinitesimal generator Q^"^' of the scaled process 



IS given as 



k=0 

oo h ^ 

+ A5")n\yj;[0(x,y + ^) -0(x,?/)]^^"^(fc), (4.2) 



fc=0 

where ix,y) £ S^ x sj"^ and (p £ BM([l,oo) x M+). 

We note that a key difference between the generators Q^^' above and A^^' in (2.2) is the 
extra factor of n in the first term of (4.2), which says that, for large n, the catalyst dynamics 
are much faster than that of the reactant. 

We will show in Theorem 4.2 that the reactant population process Y'^'^' can be well ap- 
proximated by the one dimensional diffusion Y in (4.1). Once again, the result provides a 
model reduction that is potentially useful for simulations and also for a general qualitative 
understanding of reactant dynamics near criticality. 

Theorem 4.2. Suppose Conditions 2.1, 3.1, and 3.2 hold. Then, as n ^- oo, y(") converges 
weakly in Z)(M+ : M+) to the process Y . 

We will prove the above theorem in Section 8. 



5. Auxiliary Results 



In this section we collect several auxiliary results, which will be used in the proofs of our main 
results. Recall that the quadratic covariation (or bracket process) of two semimartingales ^ 
and C is the process {[C7C]t}t6K+ defined by 

t 



K, C]t ■■= ^tCt - [ ^s-dCs - [ Cs-d^s, t > 
Jo Jo 



where ^o- '■= 0, ("o- '■= 0. The predictable quadratic covariation of ^ and C is the unique 
predictable process {(^, C)t}teiR+ such that {[^, C]t — (C; C)t}teM+ is a local martingale. If ^ = C, 
then [^] = [^, ^] and (^) = (^, ^) are, respectively, the quadratic and predictable quadratic 
variation processes of ^. 
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For X = (xi, X2, X3) G W, let </>j(x) = Xj, i = 1, 2, 3, and h := (pi — ^3. Note that for a locally 
bounded measurable function / on W 

Mt\f) := /(Wi")) - /(Wj")) - riW/(wW)ds, t > 0, (5.1) 

Jo 

is a local martingale with respect to the filtration (j(Ws : s < t). For the rest of the paper, 
we suppress the filtration, and simply refer to M^"''{f) as a local martingale. 

Let 

#):=A!")n;x!")(0)j[|\^^(„,^^^d.. (5.2) 

This process will play the role of the reflection term in the dynamics of the catalyst, arising from 
the controlled immigration. The following tightness result will be used in the weak convergence 
proofs. 

Proposition 5.1. Suppose Conditions 2.1 and 3.1 hold. Then the family {{X'^'^''\Y^'^\fj^'^^)}n<m 
is tight in Z)(M_|_ : [l,oo) x M+ x ]R_|_). // additionally Condition 3.2 holds, then the family 
{(^i+U2 - r}i"^)}neN,.eM+ is tight m D{R+ : [l,oo) x M+). 

The proof of Proposition 5.1 will be based on the following results. Lemma 5.1 below gives 
some useful representations for the catalyst and reactant processes. Lemmas 5.2 - 5.5 and 
Corollary 5.1 provide moment bounds that are useful for arguing tightness. Proofs of these 
results are given in Section 5.1. 

Lemma 5.1. Suppose Condition 2.1 (i) holds. Then the process {X^'^\Y^'^>) can be represented 
as 

Xi") = 4") + cS^^aS") /* X^r^ds + Mi")(,^3) + M,(")(/i) + f^P 

Jo 

= r ("xi") + 4"Ui") f Xi")ds + M^") {<P3) + M.(") ih)\ (t) (5.3) 

and 

ft 



yi-) ^ y(n) ^ ^(n)^(n) f ^{n)y^(n)^^ + Mt\<P2). (5.4) 

Jo 



Moreover, for t >0, 



(MW(,^3))t = Ai")ai") / xi^Us, (5.5) 



t 

t 




(MW(/i))i = A;"VS"^(0) I xWl^^(„)^^^ds, (5.6) 



and 



(M(")(</.2))t = A5")4"M xWyi")ds. (5.7) 

Jo 
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Let ifW := M(")((/)3) + M^''\h) and 

JO 

Then we have the following second moment estimate. 

Lemma 5.2. Suppose Conditions 2.1 (i) and (ii) hold. Then there is a K & (0, 00) such that 
for all neN andT>0 



E (sup [{xry + {Hrr + {Nrr + ivrr) i < exp(i^^^)(x^)^ (5.9) 

and for each A; G N 

E fsup (^1^),,)'] < eMKT'k^){y^^^f, (5.10) 

w/iere 4"^ := inf{t > : X^^"^ > k}. 

In order to study properties of invariant measures of X*-"-*, it will be convenient to allow the 
initial random variable Xq to have an arbitrary distribution on S^ . When X^ has distri- 

(n) 

bution fi on Sj^ , we will denote the corresponding probability and expectation operator by P^ 
and Efj_ , respectively, li fi = 6x for some x G S)^ , we will instead write Px and Ex , respectively. 
When considering an initial condition x for X^"', x will always be in Sj^ , although this will 
frequently be suppressed in the notation. The symbols E and P (without any subscripts) will 
correspond to the initial distribution as in Condition 2.1. 

Lemma 5.3. Suppose Conditions 2.1 (i) and (ii), 3.1, and 3.2 hold. Then there exist S,p G 
(0, oo) such that for every M > 

sup Ex I sup e''-^'"' I =: d{S,p,M) < oo. (5.11) 

neN,x<Af \0<t<p J 

Lemma 5.4. Suppose Conditions 2.1 (i) and (ii), 3.1, and 3.2 hold. Then there exist 5,d £ 
(0, oo) such that for every x G Sj^ ' ''^ ^ N; o,'>T'd t > 

^x(e5^'"')<dV^ (5.12) 



The following is immediate from Lemmas 5.2 and 5.4. 

Corollary 5.1. Suppose Conditions 2.1 (i) and (ii), 3.1, and 3.2 hold. Let 5 be as in Lemma 
5.4 and T G M+. Then there exist a d{6,T) G (0,oo) such that for all x G S^ and n G N 

sup£;J sup {Xi"^)A<d{5,T)e^\ (5.13) 

\ s<u<s+T I 
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The next lemma follows on combining Lemma 5.4 with arguments as in the proof of Lemma 
5.2. The proof is omitted. 

Lemma 5.5. Suppose Conditions 2.1 (i) and (ii), 3.1, and 3.2 hold. Let 5 he as in Lemma 
5.4. Then for each T > there are Lt, Lt G (0, 00) such that for all n gN and s G M+ 



E f sup ((Xi:l - X("))2 + (Hill - H^^f + (7vi5 - N^'^^f + {tl - fif 

\t<T ^ 



J") 



<LTE(x^"-n <LT[e^''° ). (5.14) 



In order to prove weak convergence results for the scaled catalyst and reactant processes, 
we will need to argue that the limit processes are continuous, which will be a consequence 
of the following bounds on the jumps. The somewhat stronger estimate on the jumps of the 
catalyst population in (5.16), below, will be used in the stochastic averaging argument in the 
proof of Theorem 4.2. Recall that for a process {,^f }tgK+ the jump at instant t > is defined 
as A^i := e* - it- and A^o := 0. 

Lemma 5.6. Suppose Condition 2.1 holds. Fix T, e > 0. Then, as n ^ 00, 

P ( sup (\Axi''^\ + |Ay/")|) > e ) ^ 0. (5.15) 

\0<t<T ^ ^ J 

If additionally Conditions 3.1 and 3.2 hold, then, as n ^- 00, 

sup P I sup lAX^+^J > e I ^0. (5.16) 

seK+ \ 0<t<T / 



5.1. Proofs of Auxiliary Results 

In this section we prove the results stated in Section 5. We begin with the proofs of Lemmas 
5.1 - 5.5. Using these results, we will then prove Proposition 5.1. The proof of Lemma 5.6 is 
given at the end. 

Proof of Lemma 5.1. Recall that W(") = (X("),y("), Z(")) and that for x = (xi,X2,X3) G W, 
0j(x) = Xi, i = 1, 2, 3, and h := (pi — (j)^. From (5.1) 

Zf) = 03(W|")) = 4") + /*i(")03(w("))ds + m(")(<A3). (5.17) 

Jo 

Using (2.2), we get 

00 
iW03(wl")) = Ai^^nX^") J](A; - l)^i"^(A:) = ci^^A^^^xf ^ (5.18) 

A:=0 
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and 



iW0i(wl")) = Ai")n2l(")J] 



A:=0 



Zin)^k-l 



n 



Zl 



(n) 



l^t\k) 



Ai^^n^xf ) Y. 



k=0 



k-l 



n 



2zf^ + 



k-\ 



n 



a4"^(A:) = 2\^^^xt^zt^S^'^ + Ai"^x("^a["\ 



A routine calculation then shows (see [8], Lemma 3.1.3) that 

(MW(</)3))t = AS")ai") ( Xf^ds. 

Jo 



(5.19) 



Next, since Xq = Z^ , we have 



J^{n) _ ^(n) ^ /j(wf')) = C A^"'^hi^f^)ds + Mi")(/l), 

h 
and, once more using (2.2), 

i(«)/i(w) = Ai")n/.;")(0)l|,=i}, w = (x,y,z). 

Thus with f]^^> as in (5.2), we get 

Let [/(") := 1(") - Z("). Then 

2 /■* 



(5.20) 



U, 



(n) 



2 / ut^dUi''^ + [?/(") ]i = 2 /" Ui'^^dfii'''> + 2 I U^J']dM^'^\h) + [?/(")] 
JO Jo Jo 



Jo 
where L'"' is a local martingale. On the other hand, using the generator AS""', we find 

([/^("))' = /i2(w(")) = f' M^)h\W^:'^)ds + M,(")(/i2), 
where M^'^>{h?) is a local martingale. The last two representations of ( fT/" 1 imply that 

(;/("))*= / i(")/i2(w("))ds-2 / i/i")d?7^^ (5.21) 

JO JO 

Now, with w = {x,y,z), 

iW/i2(w) = Af)nVi"^(0) n 1 - 2 + ^ 

l{x=i}- 



fl-^)' 1 



{x=i} 



AS"VS"^(0)(i + 2(l-.)V{.=i}- 
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Thus 



i(")/i2(w("))ds = / f- + 2(1 - Zi"))") dfi^J'\ 



Moreover, recalling that ?7^"^ = X^"^ — Z^'^' and that 77^"^ increases only when X^'^' = 1, we 
have 

JO Jo 

Combining the last two equations with (5.21), we get 

Also, since fj^"^' is a continuous process with bounded variation, [fj^^'] = [?}'"', M^"^(/i)] = 0, 
and consequently from (5.20), 

{M(^\h))t = {U^% = Ai"Vf^(0)^*^i"H^^,„,^^^ds. (5.22) 

Next, by (5.17), (5.18), and (5.20), we have 

X,(") = 4") + cf) aS") /* xi^Us + Mi") (,^3) + M,(") (/i) + r)!") . (5.23) 



Since n^") is non-decreasing and fi^ ^riy(n),,^d7]s = 0, we have from the characterization 

^ i^S 7^1} 

given above (2.4) that 



xfUr li 



,") + cS"Ui") I" xi^Us + M.("H</.3) + M^''\h)\ (t). 



Finally, for the reactant population, using similar calculations as for X^^\ we get 



and (M('^)((/)2))t = A^")^"^ /o ^i"^^i"^ds- Details are omitted. D 

Proof of Lemma 5.2. Recall that i/(") = M'^'^\(f)^) + M^^\h) is a local martingale, and, us- 
ing (5.5), (5.6), and the Kunita-Watanabe inequality (see e.g. [15]), we have a bound for its 
predictable quadratic variation: 



{h'^% < (^(M(")(</.3))i + ^(MW(/i))i) < 2X'r\at'^ + ^^'^^\0)) / xi^Us. (5.24) 
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Combining this with Doob's inequahty, we have 

E [sup(i7fV] < 8AS")(aS") + ^,t\^))E {[ ^'i'^^s^ ■ (5-25) 

Next, from (5.3), X^ = T{N ){t). The Lipschitz continuity of the Skorohod map impUes 

sup ixi""^ - 1| < 2 sup |ivi"^ - 1|. (5.26) 

t<T t<T 

Letting |X(")|^ j, := supj<j, iXj,^"^]^, we now get 

|^^"^|^,T < d^^"^ - M*,T + if < 2|X(") - 1\It + 2 < 8|iV(") - 1\It + 2 < 16|iV(")|2_r + 18. 
Combining this with (5.8) and (5.25), we obtain 

E (|^^"^|*,t) < 18 + WE (|iV("'|^_^) (5.27) 



< 18 + 48 



i?(xJ"V+ (r(ci"^Ai"V + 8AS")(ai") +/iS")(0))) j^ E{\x'-:^\l,)ds 



Using Gronwall's inequahty, we get, since E{Xq )'^ = Xq > 1, 

i?(|xW|2^^)<66(xi"))2exp«^), 

where i^g := 48r (t (ci^^Ai"))' + 8Ai"Vf^ + l^t\^))) ■ Since ct\x^^\ and a^ converge 
as n — )• oo, we have that for some K G (0, oo) and all n G N 

E(sup(l("))2) < exp(i^r2)(x("V. (5.28) 

t<T 

Using (5.28) in (5.25), (5.27), and (5.23), we have the estimate in (5.9) by choosing K suffi- 
ciently large. 

We next establish (5.10). Using Doob's inequality once more, we have 

Combining this with (5.7), we obtain 

E Uv (^|),,('^2))'') < 4A(")a(")i? (f^' "^xWyWd.") . 
Using (5.4), we now get 

E{\Y^X,T..[-^) ^ 3 ((^r)^ + [T{c^^h'?kf + AX^-^a^-h] £ E{\Y^-Xw^^)^') ' 

The estimate in (5.10) now follows by choosing K sufficiently large and applying Gronwall's 
inequality. D 
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Proof of Lemma 5.3. Using Conditions 2.1 (i) and (ii), 3.1, and 3.2, and using a Taylor series 
expansion, we see that there are (5o,rfi,rf2 G (0,oo) such that for aU 5 € [0, Jq] and n G N 



oo 

-5d2 <^n^ \e^^^ - 1 /iS"^(A:) < -5d 



k=0 



For 6o as above and S < 5o, let 



.(") ._ ^J 



X^)i 



and 



Note that, for any t > 0, 



fc=i 



t oo 






(fc-l)i5 

e 1 — 1 



k=0 



,W/ 



^i (^)j l{x("'>i}^^- 



and 



< af^ < e^5. 



We now argue that 

is a local martingale. Let /(x) = e and 



M,^")''' := exp (<5xf ) - /Sf )'^) - aj^^?;") 



/(x) 
Here C^^' is the generator of X''"\ that is, for x E §^ 



£W/(x)=aS"Vx^ 



fc=0 



/|(^ + ^^-^-l) +1)- /(^) 



n 



/.("^(A:). 



Note that 



and thus 



oo 

,(x) = n^A(")xj:([e^-l];.;")(fc))l 



{x>l}, 



k=0 



Kxi"))d. = /3j")'^ 



Also, 



(5.29) 



-Sd2X^r^ I Xi")l^^(„)^^^d5 < /3j")'' < -6d,\P I X^,-h^^,^„,^^^ds (5.30) 



(5.31) 



(5.32) 



(5.33) 



(n)An) 



£(")/(Xi"))l^^(„)^^^d. = af^r}} 



(5.34) 
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Consider the Markov process V^^' defined by 

V^^ := (x("\exp [- j\{n-^)dsy^ , t > 0. 

Denote by C^"^' the generator of V^"^' . Then the action of the generator on the function f{x)g{y) 
with /(x) = e and g{y) = y is given by 

C^''Hf{x)g{y)) = £(")/(x) - q{x)f{x) = £W/(x)l|,=i|. 

Consequently, 



jx 



v>(") 



exp - ; q 



-(Xi"))d.) - ^*£W/(xW)l^^(„)^^^ds, t > 0, 



is a local martingale. Using (5.33) and (5.34) we now see that the last expression equals M^ , 
t > 0, which is thus a local martingale. 

We next show that for every M > 0, S < 6o, and t > 

d3{6,t,M) := sup sup^a; (e^^^" ) < cx). (5.35) 

x<Mn&N ^ ' 

From (5.30), pf' < 0, and consequently 

yt^'' + af r)(")) e^'"" < Ali^^'' + a(")#\ (5.36) 

Thus, using (5.31) and that M^"-''^ is a local martingale, we get 



The estimate in (5.35) now follows on combining the above inequality with Lemma 5.2. 
Fix M > 0, X < M, and (5 < ^. Then, from (5.36), for all p > 0, 

^a. sup e°^« < ^x- sup Mf' +a^ 'i]l 
\0<i<p / Vo<t<p ^ 



< E^ sup Mf^-'' + e^5E^ sup r/f^ • (5.37) 

\0<t<p / \0<t<p / 



Moreover, 



0<t<p 



ii;. ( sup Mf )'^ < 4ii;,. ( (m(")'^)2) < 8 



'E^(e^^^'^'--^'^'') + {e'5fE.{fj^;)) 
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and, from (5.30) and (5.35), 



EAe 



,254"^-2/3(")'^l<,^ , . 



45xi"'^^2 ( ^ /^^_4/3<"''^\\2 



EAe 



(^ 



iW. 



(n) 



< {d3{45, p,M))^E,\ exp 4^2 A^V sup ^t j j ■ 

Choose p < ( 8d2 sup„gf^ X^ j . Then, by combining the above estimates, we can find a 
d4{6, p, M) < oo such that for ah x < M, n£N and 5 < ^ 



Ex sup e 

\0<t<p 



^^'"^ I <d4{5,p,M)Ex f exp (45^2 aS"V sup Xj^"^ 



0<t<p 



<d4{5,p,M)Ex (expf- sup xj"^ ) 1 <d4i5,p,M) 



Dividing both sides by Ex i supo<«p e 

for any x < M and n G N. The result follows. 



2 o<t<p 



<y(n) 



Ex sup e 

\0<t<p 



<5X("' 



yields 



^x ( supo<i<p e 



<5xi"' 



<(i4((5,p,M), 
D 



Proof of Lemma 5.4- For 5 G (0, 1), n G N, define 

oo 

6(")'\x):=Ai"Vxj;(e^^-l)M")(A:), 

fc=0 

oo 

6(")'2(a;) := Ai")n\ j;^ (e-^V _ l) ^^(A;), and 
fc=i 
5i")(x) := 6p(x)l^,>i| + 6f '^x)l^,=ij. 

From (5.29), we have, for some k G (0,oo), 

supbP'\x) < -6dix inf A(") < -6kx < -6k 



n6N 



n6N 



for all (5 < 5o (with 6o as above (5.29)) and x > 1. Observing that with /(re) = e , — .A 
b^ (x), where £^"^ is the generator of X^^' defined in (5.32), we have that 



uir^)._^^5xr-fobrm">)ds ^>o 



(5.38) 



is a local martingale. Fix 5 and p as in the statement of Lemma 5.3. Without loss of generality, 
we can assume that 6 < 5q. Note that on the set 

{u : xW(u;) > 1 for all s G [(j - l)p,jp)}, 
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we have 

^[^(n) _ ^(n)^^^j < ^^^(„) _ ^(n)^^^j _ / bf)^xin)^ds - 5.p = vf - 5.p. (5.39) 

•^(i-i)p 
Fix t > and let A^ G N be such that [N — l)p < t < N p. Then, similarly, on the set 

{uj : Xi"^(w) > 1 for all s G [(iV - l)p,i)}, 

5[xf)-X{;)_,);<.(;)(t), where 

"■J >P 7(j_i)p 

Now, for a fixed cj, let m = m(uj) be such that [{m — l)p,mp) is the last interval in which 
Xi"-> visits 1 before time Np. We set m = if 1 is not visited before time Np. We distinguish 
between the cases 0<m<N,m = N, and m = 0, where the latter corresponds to the case 
where 1 is not visited before time Np. 
Case 1: < m < iV. 
In this case 

N-l 



For j G N, let 
and 



6xi-^ < SX^4 + Y: i^f^ - 5-P) + 4\t). 

j=m+l 

7f^=inf{t>(j-l)Hx(") = l}Aj> 



}f:= sup[X;;) „ -X^t^■ (5-40) 

■' 0<t<p i^+Tj )^3P "/j 



" (n) (n) 

Then SXmp < S9m + <5- Combining the above estimates, we have 

N~l 

5xt^ < 56^^^ +5+ Y. i^f - ^^P) + ^N(t)- (5-41) 

j=m+l 

Thus, in this case 



,5X/") < <5X("^ + max ( 



N~l 

Y {vf^ - 5np) + 5e\-^ 
j=i+i 



}+vi:\t), 



where by convention J2i=i+ii''^i ~ ^^P) = iov I = N — 1, N , and 9q := 0. 

Case 2: m = 0. 

In this case, 1 is not visited before time Np and thus 

Af-l 

dlf) < 5xt^ + Y ^^f - ^^P) + ^S^ it) 

< 6xlr^ + max^ | ^ ^f " ^^P) + ^^f"^ | + -Pi^). 
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Case 3: m = N. 



Suppose first that there is an s G [{N — 1)p,t\ such that Xg = 1. It then follows that 



^W 



5X1''^ < 56^ + 5. 



Now suppose that there is no such s G [{N — l)p, t\. Define m' € {1, 2, . . . , A^ — 1} to be such 
that [{m' — l)p, m' p) is the last interval in which X^^' visits 1 before {N — l)p. Once again we 
set m' = if there is no such interval. 



If in' = 0, we get exactly as in Case 2 that 

N-l 



5Xr^ < SXlr' + max ^ (.f ^ - S.p) + 



Jn) 



0<1<N 



J=l+^ 



}+vi;\t). 



li 1< m' < N - 1, then 



N-l 



sxi'^use':)+s+ y: {vr-5np)+v^-\t) 

j=m'+l 

N-l 



< 6X1,"^ + max ^ {vf^ - Snp) + 59^ 



Jn) 



0<1<N 



,i='+i 



}+vi;\t). 



Combining the three cases, we have 



6X^'^^ < max <J 6XJ'''> + max { 



Thus, for any Mq > 0, 



N-l 

Y^{vf-5Kp)+5e\'''^ 
j=l+l 



^+„W(t),5 + 50H 



N-l I N-l 

< E ^- ^N\t) + E (^f ^ - ^^P) + ^^'(^'^ + ^ + ^^^^ >Mo\+ P.i^oP + 6 > Mo) 
1=0 \ j=i+i 



< g<5(i+x)-Afo y^ 



N-l 



1=0 



Ex i exp 



7V-1 

j=i+i 



-5K.p{N-l-l) 



+ Ex{e<')]. 



Recalling U^^' from (5.38) and using its martingale property, we get 



7V-1 



< e-^°e(i+^)^d(5, p, 1) (l + YZ^) ' (5-42) 
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where the last inequahty follows from Lemma 5.3 and the observation that 



g(") N . „ / ^x'"-* 



snvE^{e^"i )<sup£;i sup e"^^* < (i(5, p, 1) < oo, (5.43) 

nSN n<m \0<t<p j 

where 0^ is as in (5.40). Finally, from (5.42), we get that for all t > and n G N 
E, (ei^'"^) = y" P, (jxf ) > 21n(y)) dy 

< 1 + e(^+^)^(i((5,/>, 1) (l + _ ^^^ \ r e-2Wy)rfy < Je^^^ (5.44) 

where d=l + e^d{5, p, 1) (l + yz^W) • The result follows. D 



Proof of Proposition 5.1. We will first consider the second part of the proposition. We begin 
by showing that {Ns+- ~ ^s ss,n is tight. For that, in view of (5.14), it suffices to show that 
the following condition (Aldous-Kurtz criterion) holds: For each M > 0, e > 0, and 7 > there 
are ^o > and no such that for all stopping times {T„}„gpj with r„ < M, we have 

sup sup P(|iv(:;t„+e - ^i+l I > 7) < e. (5.45) 

,n>ni) 9<So 



Let M, e,7 E (0, 00) be given. Note that 

CS+T„+9 
s+r„ 

By (5.13) we have, for 60 sufficiently small, 

rs+T„+d 
seIR+,neN6l<(5o 

It remains to prove that, for some 5q > 0, 



^(l^il+. - ^ill > 7) < ^(|cS"^aS") r^^'xHdnl > |) + P(|<1„^, - <lj > |). 
r (5o sufficiently small, 
sup sup P f 14") Af) Z''""^' Xt^du\ > j) < ^. 



sup sup P (|<l„+, - H^:l^ >?)<!• (5-46) 

Using the martingale property of H^"^' , 
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and, using the properties of the predictable quadratic variation, 

(, , , , fS+Tn+B , s / s rS + Tn+d 



Now, using (5.13) once more, we can choose Jq > such that (5.46) holds. This proves tightness 
of {Nil. ~ ^s }s,n and, using the continuity property of the Skorohod map (from L>i(R+ : M) 
to DiR+ : [l,oo))), that of {xif - li"^},,„ and {fiif - r?i"^},,„. Tightness of {xil\}s,n now 
follows by using the uniform estimate in Lemma 5.4. 

Now we consider the first part of the proposition. Tightness of {X^^\fj''"') follows as before. 
We now consider Y^"^' . Fix e > 0. Using (5.10), we get, for K G (0,oo), 



P(sup(y/")) >k)< P(sup(y(("„)^p > K and af > t) + P (4"^ < t) 

P(sup,<r(i;tK/) ^ ^^^^^ ^ P.nr7<'T2i.2v>)^2 Krsun..^rX^)^2^ 



+ p(sup(xH) > k) < -P(^^:^^)(^r)- ^ Eisu,,,AXl 



Using (5.9), we can choose k such that 

P(sup,<r(xf))2) ^ 

sup^ T2 -<2- (^-^^^ 

Now choose K such that 

exp(Kr2yfc2)(y("))2 g 
sup ^- ^^ < - . 

n6N K^ 2 

The last two displays imply sup„gt^P(supj<2-(y/" ) > K) < e, and since e > is arbitrary, the 
tightness of the random variables {Y^ }n£fq, for each t > 0, follows. To establish the tightness 
of the processes {y^^jneN, it now suffices to show that for each M > 0, e > 0, and 7 > there 
are 6q > and no such that for all stopping times {T„}„gN with r„ < M, we have 

sup sup P (itite - i;L"^| > 7) < e- (5.48) 

n>no 0<So ^ ' 

Fix M,e,7 E (0,oo). Then, for any Q G (0, 1), 

p (lei. - «:'i >- 7) < ^(ic...,i-' - 1:« I ^ ^) + ^ ^ s '' + ') ■ 

Taking T = M + 1 and /c as in (5.47), we have P{(yf' < M + 1) < e/2 for all n G N. 
For the first term on the right hand side of the last display, we get, using (5.4) and that 
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S^Pt<TAal"' ^t"^ ^ ^' 



p I ly^"'' 






Jr^Aa^") 2A;/ V (r„+e)A4 TnAa)^ 2 



(„,(<^2)-M("J („)(<^2)|> 



(5.49) 

7^ 



The first term on the right hand side can be bounded as follows: 



•'2 ^2 



(T„+9)Aa^"' 

. (n) 
TnAcr^ 



yi")df 



2^7 ~ V 7 



^(") 



>i^<r?^££^VErrr"""~""'«»).r 



T,iA(T^ 



< 



,(") \ (") \ 2 



AKCr\ A 



2 ^'2 



7 



exp {K{M + l)2fe2) (y^' 



2j,2\ /,,Ws2 



where i^ is the constant from (5.10). Thus, for 5q sufficiently small, we get 



sup sup P 

neN6/<(5o 



JTnhcrl. 



^i H'/^- 



The second term on the right hand side of (5.49) can be bounded as follows: 



P(\M 



in) 



(t„+6»)A(t^' 



,„,(..)-M;-i,,,,.,|.2). "(<"''''^'>>'- -^->-r'-<-'"''..)W, 



Tnha), 



(7/2)^ 



{T„+6l)Aa^' 



(") 



TnAo-J, 



(n) 



xi")yi")d5 1 < 4A5"^a^"^A:ei5; ( sup y("\„, ) 

y ~ 7^ \o<.<A/+i «A4")y 



Using (5.10) once more, we have that, for 6q sufficiently small, the second term in (5.49) is 
bounded by |. Combining the above estimates, we now see that (5.48) holds, and thus tightness 
of {y(")},,6N follows. D 

Proof of Lemma 5.6. Consider (5.16). Let N^j- be the number of deaths of particles of the 
(unsealed) process X^"^' in the time interval [s, s + T]. Fix e, (5 > 0. Then 



.(n)| 



.(n)| 



.(n) 



.(n) 



P( sup |AX]:;;| > e) < P( sup |AX]^;| > e; sup X'^;^^ < nL) + P[ sup X"^;^, > nL). 

0<t<T 0<t<T 0<t<T 0<t<T 



By Corollary 5.1 we can choose L G (0, oo) such that P(supo<t<T-'^s+t > ''^-^) < 3, for s G 
and n G N. Next, consider 



^(")| 



.(n) 



P{ sup |AX,7,| > e; sup X^^ < nL) 

0<t<T 0<t<T 



.{n)| 



rW 



.(n) 



rW 



< P{ sup |AX^"^J > e;N^'tl < nCL) + P{ sup X,7j < nL-N'^J'^ > nCL). 

0<t<T ' 0<t<T 
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Note that on the set {supo<t<T X^t — '^-^1 the branching rates of X'"-* are bounded during 
the time interval [s, s + T], uniformly in s and n, and thus we can choose a C G (0, oo) such 
that for s G M+ and n G N 

P( sup Xi:;i < nL;N^''^ > nCL) < -. 
0<t<T ' 3 

Finally, let, for n G N, {Q" lieN be i.i.d. random variables distributed as fJ-i ■ Then, since the 
variance of the offspring distribution converges, we have for uq sufficiently large and all n > uq, 

P \ sup I a1,^"^,| > e; N^Jll < nCL ] <P\ max *^' ~ ^* > e ) 

\0<t<T ' J \l<i<nCL n J 



. E^i^r^-ii^-)^ g ^;.,. ^ <i 



j=l i=l 



Combining the above estimates, (5.16) follows. The limit in (5.15) can be established similarly, 
using Lemma 5.2 instead of Corollary 5.1; the proof is therefore omitted. D 



6. Proof of Theorem 2.1 



The following martingale characterization result will be useful in the proof of Theorem 2.1. 
The proof is standard and is omitted ([17, 11, 12], see also Theorem 5.3 of [2]). 

For (p G C^°°([l,oo) X M+), let 

d I d^ d 1 d^ 

C<j){x, y) := ciXix—(j){x, y) + -ai\ix-—^(j){x, y) + C2\2xy-^4'{x, y) + -a2\2xy-g-^4>{x, y). 

Let Cl := Z?(]R_|_ : [l,oo) x M^) and J- be the corresponding Borel c-field (with respect to the 
Skorohod topology). Denote by {Ft}t&,+ the canonical filtration on {(l,F), i.e. Ft = (T{Tis\s < 
t), where 7rs(a;) = Cos = oJ{s) for cD G f^. Finally, let 7r(*\ i = 1,2, 3, be the coordinate processes, 

i.e. {TTW{u),TT(^\io),7r^^\u;)) = TT{u). 

Theorem 6.1. Let P be a probability measure on {^,F) under which the following hold a.s.: 

(i) TT^^' is a non- decreasing, continuous process, and ttq = 0, 
(a) [tt^^' ,7r^'^') is an ([l,oo) x M4.) valued continuous process, 

(Hi) /o°° l(i^oo)(7ri^^)d7rf ^ = 0, 
(iv) for all (f) G C^([l,oo) x M+) 

is an {Ft} martingale, and 
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(v) P o (ttq , ttq )~^ = P o (Xq, Yq)~^ , where X, Y, and P are as in Proposition 2.1. 
T/ien Po(7r(i),7r(2))-i = Po{X,YY^. 
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Proof of Theorem 2.1. Recall that for cj) G C^{[1, oo) x M+), we have 



where M^ {(f)) is a martingale, and A}-"^' is as defined in (2.2). Also note that A^"^' can be 
rewritten as 

iW</.(x,2/) = £(")<^(x,y) +pW<A(y)nAS"VS"^(0)l{.=i 



where 



C^^)4){x,y):=Xf^n'xY, 



k=0 

oo 

+ Ag 'T. xy >^ 
fe=o 



x-\ ,y ] - (t){x,y) 



n 



{x=l}, 



,W 



f^rik) 



x,y + 



k-1 



n 



4'ix,y) 



l^t\k) 



and 



pW</.(y):=n(^<A(l,y)-,A(l-^,y' 
Thus, using (5.2), (6.1) can be rewritten as 

JO 

+ rpW</.(yW)d,7W + Mi")(<^). 



Recall the path space (i^,J^) introduced above Theorem 6.1. Denote by P'") the measure 
induced by {X^'^\Y^'^\'r)^'^') on {(l,F) and by £"("■) the corresponding expectation. 

From Proposition 5.1, P'"^ is tight. Let P be a limit point of {P*^"^} along some subsequence 
{rifc}. In order to complete the proof, it suffices to show that under P properties (i) - (v) in 
Theorem 6.1 hold almost surely. Property (i) is immediate from the fact that 1)^'^' is non- 
decreasing and continuous with initial value for each n. Also, property (v) is immediate from 
the fact that {Xq ,Yq ) = (1, 1), a.s., for each n. Next, consider property (ii). The continuity 
of 7r(i) and vr^^) follows by (5.15) (see [7], Proposition VI.3.26, p. 315). 



To see (iii), consider, for (5 > 0, continuous bounded test functions /^ : [l,oo) 
that 



fs{x) 



1, ifx> 1 + 2(5 
0, if X < 1 + 5. 



such 



(6.2) 
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Note that, for each n G N, /q°° fs{Xs )dfjs = and thus, for each 5 > 0, 

Consequently, for each 6 > 0, Jg°° l[i+2<5,oo)(^s )dTTs = 0, almost surely w.r.t. P. The property 
in (hi) now follows on sending (5 — )• 0. 

Finally, we consider part (iv). It suffices to show that for every < s < t < oo 

where -0 : fi — t- M is an arbitrary bounded, continuous, Tg measurable map. Now fix such s,i, 
and ip. Then by weak convergence of P^^*'i to P and using the moment bound in Lemma 5.2, 

hm ^("^■)(V^(.)(0(vr«,4^V<^(-F,vrf )- f C^P M^^)du - T ^(1, vr(^))d4'))) 



fc— >oo 



dx 



^(^(.)(<A(vr«,vrf))-<^(vrW,vrf))- / CHn^^\n(J^)du - f ^(1, ^^^M^^))). 



dx 

To complete the proof, it suffices to show that the limit on the left side above is 0. In view of 
the martingale property in (6.1), to show this, it suffices to prove that for (j) E C^([l, c«) x M_|_) 



and 



lim E 

n— >oo 



lim E 

n— >oo 



£(n)0(xW,yW) _ /:^(xi"),i;("))) ds 







j^*(z^w^(yW)-|^(i,yW))dr)W 



(6.3) 



(6.4) 



The latter is immediate upon using the smoothness of ip and the moment estimate for i)^'^' in 
(5.9). For (6.3), we rewrite C^'^'cp using a Taylor expansion as follows: 



/:W0(x,y) = AS"^n\^ 



fc=0 



+ X2 iT-^xyy^^ 

k=0 

dx 



k-l d 
n dx 

k-l d 



1 n-iV 52 



n dy 



, , 1 /fc-iV 92 ^. ■ 

(a;, y) + o ( — ;— ) -^^^\^^ y) 



2 \ n J dy 



f,^^\k)+R^-\x,y) 



c\'''Xl'''x^(p{x,y) + -apApx— ^(/>(x,y) 



+ cth'^^xy^^{x,y) + ia^^h^^^xy-^cP{x,y) + R^^\x,y), 



»xW.„,_^(^^y) + 2"2"^^2"^'^y^' 



where the term i?'"'(x,y) is a remainder term, which, using part (iii) of Condition 2.1, is seen 
to satisfy sup|2,| |y|</^ |i?("'^(x,y)| — ;• as n — )• 00, for any L G (0,oo). Furthermore, using the 
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compact support property of </), it follows that Imin^oo E j^ \w-^''\Xs ,Ys )\ds = 0. Next 
note that 



£("V(x,y)-ii(")(x,y)-£0(x,y) 



(Apcp - \ici)x—(l){x,y) + -(Apap - \iai)x-—^(i){x,y) 



+ (4"^ 4"^ - >^2C2)xy-Q-(l){x,y) + -(Ap4 - X2a2)xy-^(j)ix,y), 



and therefore, in view of Condition 2.1, 



sup 



C^''Uix,y)-R^''\x,y)-C^ix,y) 



as n — )■ oo. 



Once more using the compact support property of (f>, it follows that 
Combining the above estimates, we have (6.3), and the result follows. 



ds = 0. 



n 



7. Proofs of Results from Section 3 



7. 1 . Proof of Proposition 3. 1 



Uniqueness of the invariant measure of X is an immediate consequence of the non-degeneracy 
of the diffusion coefficient (note that a2\2X > a2A2 > 0). For existence, we will apply an 
extension of the well known Echeverria criterion for invariant measures of Markov processes 
([12], see also Theorem 5.7 of [2]). This criterion, in the current context, says that in order to 
establish that a probability measure ui is an invariant measure for X, it suffices to verify that 
for some C > and ah cp G C^°°([l,oo)) 



Ci(j){x)i'i{dx) + CaiAi(/)'(l) = 0, 



where 



l,oo) 



Ci<j){x) = ciXiX(j)'{x) H — aiXiX(j)"{x). 



(7.1) 



(7.2) 



We now show that (7.1) holds with vi = vi and C = ^%^. For (f) G C^([l, oo)) and p as in 
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(3.1), 

ci\ix(f) {x) + -aiXixcp" {x) I p{x)dx 
1 \ ^ / 

oo 



ciXiOe '^^^(pix] 



2ciXi0 — e "i^^(j){x)dx 
1 Ji "1 



+ -aiXiOe °i (x) - / Q1A16' — e °i cp (x)dx 
2 1 Ji ai 



Thus (7.1) follows. D 



7.2. Proof of Theorem 3.1 



Throughout this section we assume that Conditions 2.1, 3.1, and 3.2 hold. This will not be 
explicitly noted in the statements of the results. 

Existence of a stationary distribution z^ J" oftheSj^ = {-\l G {n,n+l, . . .}} valued Markov 

process X^"'' follows from the tightness of {X^ }t>o, which is a consequence of Lemma 5.4. 
The uniqueness of the stationary distribution follows from the irreducibility of X'"'. 

In order to establish the tightness of the sequence {I'l }neN, we will use the following 
uniform in n moment stability estimate for X^""' . 

Theorem 7.1. There is a to £ M+ such that for all t > tQ and p > 

limsup4i?.(te))')=0. (7.3) 

Proof. Fix an L > 1, and let r^") := inf{t : X^^"^ < L}. Observe that if t e [{N - l)p, Np) for 
some N gN, then, following arguments as in the proof of Lemma 5.4, for x > L, 

^x(r(") >t)<pAY^ {vf^ - 5kp) >6{L-x)] < S^-L-&kp{N-i)) _ 

Thus we have that 

supP,(t(") > t) < 7ie^^e-^2*, 

nSN 

where 7^ S (0, 00), i = 1,2. The above estimate along with Lemma 5.4 implies, for n S N, 

S <>(") / S ■Cr{ri)\ / S <>(") 






<dl(l + 6-5^6-^*), 
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where d is as in Lemma 5.4 and di G (0, oo) is some constant, independent of n. Fix p > 0. 
Then, for some d2 G (0,oo), we have 

sup — < sup < sup ■ 



neN XP n&N X^ „gN xP 

Choose to large enough such that ^to > ^- Then for t >to 



hm sup = 0. 

^^°° nm XP 



The result follows. D 



As a consequence of Theorem 7.1, we have the following result. For 5 G (0, oo), define the 
return time to a compact set C C [l,oo) by t^^ (6) := inf{t > S\Xf G C}. 

Theorem 7.2. There are c,6 £ (0, oo) and a compact set C C [l,oo) such that 



fr^c\h 



sup^;^ ( / "^ ln(x["^)dt j < cx^ a; > 1. 



Proof. Note that ln(X|"' ) > since X^ > 1. Applying Theorem 7.1 with p = 3, we have that 
there is an L G (1, oo) such that with C := {x £ R_|_|x < L}, for all x G C^, 



^npE4{X^f)<lx^ (7.4) 



where to is as in Theorem 7.1. Let 5 := IqL and r^") := 4"^(<5) = inf{t > 5\xf''^ < L}. 
Consider a sequence of stopping times defined as follows: 

a("):=0, a(-):=a^:l,+to(x^t, V l) , m G N. 
Let TOq = min{m > 1|X "^j < L}, and 

V^^\x):=Ejr ln(X,("))dt J . (7.5) 

Then 
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Let Jj" := a{Xs |0 < s < t}. We claim that there is a cq G (0, cxd) such that for all n, /c G N, 
x>l 

r^"' \ / \ 3 



E^ ( / ^^^ HXDdAJ^^, l^,„^, < CO U-, 1^^^^,^. (7.6) 



fc 



Due to the strong Markov property, to prove the claim it suffices to show that for some 
Co € (0, oo) and for all n G N, x > 1 

r ln(xf))(itj <cox\ 

Note that for x > 1, o"}"" = to{x V L) < cqx, where cq = t^L. Using this bound along with 
Lemma 5.2, we get, for some cq G (0, oo), 

E, I sup ln(i;"M < In I E, I sup xiA j < In (x^e^^^"^)') < cox^. 

The claim follows. 

From the estimate (7.6), we now have 






(7.7) 



supy(")(x)<cosupSJ Y, [X%) . 

Note that \ X,L\ is a Markov chain with transition probability kernel 

P(")(x,^) := pI^I^^^{x,A), X G [l,oo), A G e([l,oo)), 

where Pj is the transition probability kernel for X^"^'. Using (7.4) and Lemma 5.4, we get 
that for any L G (1, oo) there exists a 6 G (0, oo) such that for all x G [1, oo) 



oo /•oo 



sup / y3p(n) (^^ ^y^ ^ g^p / y3p(n) ^ (^^ ^y) (73) 

/"OO /"OO -| 

= sup / 2/^P/g'2(^'^y)l{x>L} +sup / y^i^S(x,dy)l{^<L} < x^ - -x^ + 61[o,i,](x). 

n Jl n Jl ^ 

The above inequality along with Theorem XA^flfl of [14] yields 

(n) , 



/ '"0 -^ / \ 3\ / / "'0 -"- / 

= 2(x3 + 61[o,L](x)) <cx3, 



" (n) (n) 

where the equality in the last display follows from the fact that ^ („) > L for 1 < fc < rrv^ . 
The result follows on combining the last estimate with (7.7). D 
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The following theorem is proved exactly as Proposition 5.4 of [3]. The proof is omitted. 

Theorem 7.3. Let f : [l,oo) — )■ M+ be a measurable function. Define for 6 £ (0,oo) and a 
compact set C C [l,oo) 

V^^\x):=eJJ^ f{xi''^)dt\, xE[l,oo). 

If sup^gpj V^"^ is everywhere finite and uniformly bounded on C , then there exists a k £ {0, oo) 
such that for all n £N,t > 0,x £ [1, oo) 

ii?. (v^-Hxt^)) + i 1^ E, (/(xW)) ds < iy W(x) + k. 



We now return to the proof of Theorem 3.1 and establish the tightness of {vi }n&n- We will 
apply Theorem 7.3 with f{x) := ln(x), and 6,C as in Theorem 7.2. Since v^ is an invariant 
measm'e for X^"', we have for non-negative, real valued, measurable functions <& on [l,oo) 

E,(<^{xt'^))u'C\dx) = j $(x)z.i")(dx). (7.9) 

Fix A; G N and let V^^^^lx) := V^'^Xx) A k. Let 

By (7.9), we have that J^ ^^^\x)v^{'\dx) = 0. Let 

$W(x) := -y(")(x) - -E^ ('T/(")(lf M . 

By the monotone convergence theorem \I'^"' (x) — )• ^'"^(rr) as A; — )• oo. We next show that 
^^" (rr) is bounded from below for all x G [1, oo): If V^'^'{x) < k 

<\-) = ]vt\x) - \e^ {vt\xt'^)) > -v^-\x) - \e^ (FW(xf); 

where the last inequality follows from Theorem 7.3. liV^'^'''{x) > k 

^t\x) = \k-\E^{vi-\x'r^))>0. 
Thus ^^^\x) > —k for all x > 1. By Fatou's lemma, we have 

/"OO /"OO 

/ ^(") {x)vP {dx) < lim inf / ^^f^> (x)z^("^ (dx) = 0. 



> — K, 
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By Theorem 7.3 we have '^^'^'{x) > j J^ Ex if{Xs)j ds — k. Combining this with the last 
display, we have 



0> 



^(")(x)zyi"^(dx) > 



t fOO 



Jl 



EAf{xi^'^)]vt\dx)ds-k. 



Using the invariance property of i^j" once more, we see that the first term on the right side 

above equals f f{x)u\^ (dx), and therefore ^ f{x)v\^ {dx) < k. This completes the proof of 
tightness. 



,W 



Wi 



The tightness of {v\ }neN implies that every subsequence of {v\ } has a convergent sub- 
sequence. Call such a limit u^. Theorem 2.1 and the stationarity of Ui imply that u^ is a 
stationary distribution of X. Since the stationary distribution of X is unique, we have v\ = J^i, 
which completes the proof. D 



8. Proofs of Theorems 4.1 and 4.2 



8.1. Proof of Theorem 4-1 



In order to prove the result, we will verify that the assumptions of Theorem II. 1 (more precisely, 
those in the remark following Theorem II. 1) in [16], pp. 78-79, hold. For this, it suffices to show 
that for all A; e N, $ G BM(M^), 4> G C~(M+), and < ti < ^2 < • • • < tfe+i < T < oo, there 
exists a sequence /i„ with limn_>.oo hn = and 



sup 

ie[tfc+i,T] 



E 



$(!>("),, 






v(") 

^ t+h„ 



^t 



KLct> ( y/") 



o{hn), (8.1) 



where C is given as 



£(/)(y) := C2X2mxy(t>' {y) + -a2X2mxy(l)"{y), </» G C; 



Ti 



?.2) 



- (n) 

Letting X° := X., , i > 0, we see, using scaling properties of the Skorohod map and 
straight forward martingale characterization results, that X° has the same probability law as 
the process X that was introduced in Proposition 2.1 with initial value Xq = xq. The following 
uniform moment bound will be used in the proof of Theorem 4.1. 

Lemma 8.1. There exists a (5o G (0, oo), such that whenever X is as in Proposition 2.1 with 
initial value Xq = x, for some x G [1, oo), we have 



sup Ex fe^"^*) =: d{5Q,x) < oo. 



0<«OO 
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Proof of Lemma 8.1. We begin by establishing exponential moment estimates for the increase 
of X over time intervals of length Ip when the process is away from the boundary 1, where 
p > and / > 1. Fix p G (0,oo). Let a := ciAi, h := aiAi, and 5 £ (0,-f A 1). Note that in 
view of Condition 3.1, a < 0. Define ar := inf{i G [0, oo)| f^ Xsds > r} and pi^r '■= Ip A 0"^. 
Then 

E^ ( exp Ua I '"' Xsds + SVb f '"' y/jT^dBfj j 

= E^ ( exp ({6a + 6^b) I ''' X^ds + 6Vb f ''' y/X^dBf - 6^b I '"' X^ds) J 



^ JO 



< (e^^^a^^^^^Y (e^^^ (2,^ /"■' ^,Bf - <5^ /"■' x.i.' ' ' 



where the inequality follows on noting that pi^r ^ lo''"^ -^sd.s and 6 G (0, — |). Using the 
super martingale property of the stochastic exponential, we have that the second term on the 
right hand side of the last display is bounded by 1. Thus, sending r — t- 00, we have, with 
-e := 6a + S^b < 0, 

E^ (exp (da I " Xsds + 6Vb T ^/xldsA \ < e^Pi^^^^'^^ = e"^'". (8.3) 



Next, for x G [1, cxo), we have by application of Ito's formula that for t < p and 6 < 5 

E. (e-^^*) < e^" + ~5e^E^r]p < J^' + xCi{p, 6), (8.4) 



where Ci{p,5) G (0,oo) and the last inequality follows by an application of Gronwall's lemma 
and the Lipschitz property of the Skorohod map (see (2.4)). 

Using the above estimates, we will now establish certain uniform estimates on the tail 
probabilities of X/^p, which will lead to exponential moment estimates at these time points. 
Fix L > 1, and let 

Tj := inf{t > (j - l)p\Xt < L} A jp and Cj := Xjp - Xr^, j > 1, 

eo = 0. Fix fe G N, and let 



,k 



inf Xs<L\ 



M := max < j = 1, 

if there is an s G [0, kp] such that Xs < L and set M equal to otherwise. Let 

Vj := aXsds+ ^/bX^dBf , j > 1. 

J{i-i)p J(.i~i)p 
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Then Xkp = Xmp + J2j=m+i ^j- Letting Ci ■= e.i + I]j=i+i ^j' we have, using (8.3), 

fc / \ k 

P,{Xkp >K)< P^(Xmp + Yl ^J- > ^) ^ P- (o™xfc^' >K-L\< Y^PxiCi >K-L) 

j=M+l ^ -*- ^ i=0 



j=0 
k 

i=0 



i=0 



Next note that, from (8.4), 



^xe''^' < E^ (e^^''^^-^-K,<ip) +l = E, (e-^''^nr,<ipEx^. (e^^'")) + 1 



<^,: 



-SXr, 



,&Xr. 



e~'^-^ U<ip e'""-^ + Xrfiiip, 5)] +l = C2{p,5) + l. 



Hence, 



P^,{Xkp >K)< iC2{p,5) + l)le-i(^-^) ^e^^'^^ < {C2{p,5) + i; 



16 2 
2 



«(K-L) 



/=0 



l-e-2' 



The last estimate yields, analogously to (5.44), 



sup Ex{ei^''p) < Ce2^', for ah x G [l,oo), 
fceNo 



(8.5) 



for some C G (0,oo). Finally, letting Sq := ^, we have from (8.4) for t G {{k — l)p,kp], k > 1, 
i5;,(e^°^') = E,{Ex,,_,Je''>''^)) < ii;4e^°^(^-)p + X(,_i)^Ci(/.,<5)) < Ce^(l + ;^Ci(p,5)). 



The result follows. 



D 



Remark 8.1. Note that Lemma 8.1 and the scaling property noted above that lemma say that 
for all X G [1, oo) 

sup sup Ex ( e''"^'" ) < oo. 
neNO<t<oo ^ ^ 



We now prove Theorem 4.1 by showing (8.1). Let, for cj) G C'^il 



'+;> 



Cx(t>[y) ■= C2X2xy4>'{y) + -a2A2a;#"(y), {x,y) G [l,oo) x 



Then 



E 



;$(y,;-),...,y/;))(0(yit 

t+hn 



i>(")^ 






<i>(y,;"),...,i>(") 



^{Yt 



C^,^)ct>{Y^'''^)ds 



X 



^{Yi:\...xt^ 



t+h„ 



£w„),/,(yW) -£-(„,</>(!>/"))) ds 



XI 



.6) 



•7) 
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For the second term, we have, using Remark 8.1 and the fact that the function $ is bounded 
and (f) as well as its derivatives are continuous with bounded support, that 



sup E 

*6[tfc+i,r] 



(y/;), . . . , y/;)) ^*^'" (/:^(.)<^(yi")) - /:^(.)<^(i^/"^)) ds 



sup E 

<6[ife+i,r] 






^[Yr,---,yt 



tk 



t+h„ 



C2A2XW (yW0'(yW) - y/"V'(^/"^; 



+ la,x,xr> (vrriYn - YrriY, 



V-{n) (v{n)ill(Y<<n)\ _ ■{>(") w^"rv('")' 



\ds 



o{hn). 



Recalling the definition of X° above Lemma 8.1, the first expected value on the right hand 
side in (8.7) equals 



huE 



$ 



tn+hnn 



/ rinll Jin 



Thus 



E 



$ y 



v>W vC") 



- ti ; • • • ; 



^C«) 



i>H^ 



i>(")- 



YtT m'zj - ^(n'") - Kc^Yn 



E 



E 



■t, ' • • • ' -'tt 



"-n^ Jin 



tn+hnn 



^[Y^!"^ yi"M/iJ^ / '■"£xj<A(i^/"^)rfs-£0(y/")) 



+ 0{hn) 



$ (y/,"\ . . . , y/;)) K (c2A2y/"V'(>;^"^) + ^a2A2y/"V"(i^/"^)) 



"-n^ Jtn 



tn+hnn 



Xgds — mx 



+ o{hn) 



To complete the proof, it thus remains to show that for some sequence {/in} with limn-i.00 hn = 



E 



1 



hnn 



tn+hnn 



X°ds — mx 



tn 



E 



hnn 



tn+hnn 



Xgds — mx 



tn 



^.9) 



converges to uniformly in t € [tfc+i,r]. From the ergodicity of X and the moment estimate 
in Lemma 8.1, it follows that 



E 



1 /■*" 
tn Jo 



Xgds — mx 



0, as n — )■ 00. 



^.10) 



The above result, along with Lemma 8.2, below, implies that there is a sequence {/in} such 
that limn-5.00 hn = and the expression in (8.9) converges to uniformly in i G [tk+i,T]. This 
completes the proof. 



The proof of the following lemma is adapted from Lemma II. 9, p. 137, in [16]. 
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Lemma 8.2. Let < t^+i <T <oo. If for all t G [tk+i,T] 



E 



1 /■*" 

— / Xgds - mx 
tn Jo 



0, as n ^ CO, 



then there is a sequence {hn} such that /i„ — )• as n —t- oo, and 

1 ftn+hnn 



sup E 



i^n"^ Jtn 



Xsds — mx 



0, as n ^ oo. 



Proof. Let a{T) := sup^>^ ^\u lo ^sds — mx\ ■ Note that a{T) converges monotonically to 

as r — )• oo. For t G [tk+i,T] we have 



E 



"n^ Jtn 



tn+hnU 



E 



Xsds — mx 
tn + hnn 1 



tn+hnU ^^ -j^ />tn 

Xsds — - / Xgds — mx 

hnn tn + h^n Jq hnn tn Jq 



tn + hnn tn 3T 

< — aitn + hnn) + a(tn) < -—aitk+in) 

hnn hnn hn 

for ah n such that /i„ < T. Note that the right hand side of the last display is independent of 

t £ [tfc_|_i,r]. Choosing /i„ = y^a{tk+in), the lemma follows. D 



8.2. Proof of Theorem 4.2 



As in the proof of Theorem 4.1, it suffices to show that for ah A; G N, $ G BM(R^), (j) £ 
C^(M+), and < ti < ^2 < • • • < ife+i < T < oo, there exists a sequence /i„ with lim„_j.oo hn = 
and 



sup 

t£[tk+i,T] 



E 



^{ltt\...,lif 



^t + hr 



y/")) - hnL^ (y/") 



o{hn), 



where C is given as in (8.2). 

Let for 4> G Q°°(M+) and {x,y) G [l,oo) x M+ 



4"V(y):=Ar^n2xyJ] 



A:=0 



<t>[ y 



k-l 



n 



4>{y) 



in) 



/i^l^), 



and recall Cx from (8.6). Then 
E 



$(y,^"^...,y/;^)(<^(i;X)-<^(y/")) 



■y{n) v{n)\ ( j^(iy{n) 

k 

t+h„ 



(8.11) 



■-E 
+ E 



c^(y,;"),...,y/;) 



c%ct>{Yi-^)ds 



ci>(y,;"\...,y/;) 



t+/i„ 



4t)K^^^"0-4?")K^*^"0 



ds 
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Using Lemma 5.4, we get, as in (8.8), that the second term in the last display is o{hn) uniformly 
in t G [tfc+i,T]. Thus 



E 



^ (Yi-\. . . ,y/;)) Uy^X) - '^(^/"^) - hncm^^') 



(8.12) 



E 



^Y,^:\...xt^ 



t+h„ 



c 



(") A fvi'T') 



X 



{")' 



Yr -c^,„,cp Yr ds 



X. 



v>(") 



+ E 



^[Yi^\...,Yt^:^]h 



■y rt+hn 



n Jt 



viri) 



1>W' 



C^,r.)ct>{Yr)ds-Cct>{Yr)]+o{h 



Calculations similar to those in the proof of Theorem 2.1 show that the first term on the right 
hand side in the last display is o(/i„) uniformly in t € [tfc+i,T] (see proof of (6.3)), while the 
second term can be written as 



E 



0(n) v>(n) 



'A"^) A'tv'y'^)-^ 



1 



$ y,f ^ . . .,Y^:>\K{c,x,Yr<^'{Yr) + ^-a,x,Yr<t^"{Yr) 






h 



n Jt 



s — mx 



To show that the latter term is o{hn) uniformly in t G [t^+i, T], it suffices to show the following 
result. 



Theorem 8.1. As n ^ oo 



sup E 

t6[tfe+i,T] 



/ 1 rt+hn 



where /f := a{ixr ,Y}'"^) ■■s<t}, and E -(„){■) = E{- \ Pr). 

■'t 

In order to prove this theorem, we need the following three results. Let 5 := L'(M+ : [1, cxd) x 
]R_l_), V{S) be the space of probability measures on S, and, given a sequence {t.„} C [tfc+i,T], 
[In be a sequence of V(^S) valued random variables defined as follows. For A € B{S), 



f^n{A) 



1 f-ntn+nhn , •, , y 



nhn 



p[{x':::^..€-i.-fi'n^A\Tnds. 



^.13) 



ntn 



Let Sq :=C(M+ : [l,oo) x E+). 

Lemma 8.3. The family of V{S) valued random variables {/injneN is tight, and any weak 
limit point is a V{Sq) valued random variable. 



Let vr = (tt*-"^', Tr*-^') with vr'"^' and vr'^-' being the canonical coordinate processes on Sq. 

Lemma 8.4. Let ji be a weak limit point of {fin} given on some probability space {Qo,J-o,Po) 
Then for Pq a.e. uj £ Qq, ft{uj) satisfies the following: 

(a) //(u;)(7r(i)(t + -) G F) = fi{uj){'K^^^ e F), forallt>0, F e 13{C{R+ : [l,oo))), 

(b) TT^^' is non- decreasing and tTq = a.s. fJ-i^j), 

(c) /q°° l(i,oo)(7rn )d-Ku = a.s. fi{uj), and 
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(d) under ix{tjj), for all cj) G C^([l, oo)) 



vr; 



(1) ^(2)^ . 



is a {Qt] -martingale, where L\ is as in (7.2) and Qt ■= cr{{Trs ,tts ) '■ s < t}. 

We postpone the proofs of Lemmas 8.3 and 8.4 until after the proof of Theorem 8.1. The 
following is immediate from the above two lemmas, Proposition 3.1 and the martingale char- 
acterization of the probability law of the process in (2.5) (see Theorem 6.1). 

Corollary 8.1. Let {X, r]) be as in Proposition 2.1 with Xq ~ vi and vi given as in Proposition 
3.1. Let fiQ be the probability measure on Sq induced by (X, r/). Then fin converges weakly to 



Proof of Theorem 8.1. It suffices to show that for an arbitrary sequence {t„} C [tfc+i,T] we 
have, as n — )• oo, 

cntn+nhn 



E 
Since 



(1 rtn+h„ 
— / X(")(is - mx 



E^ 



(n) 






E 



ntn+nhr, 



E 



r(") 



1 



^r \nh / 



X^^Us-mx 



nt„ 



0. 



X^J'Us 



'-n Jnt 



it suffices to show that 



E 



vr^ 'diXn - T^o dfiQ 



(1)^ 
TTq a fin, 



as n — )■ oo. 



.14) 



For any c > 0, let ipc be the following continuous function: 



1, if x < § 
0, if X > c 



and tpc is linearly interpolated on [§,c]. By Corollary 8.1 fin converges weakly to fiQ, and 
therefore, for every c > 0, 



E 



■n-Q \jci-^o )dfin - / vr^^Vc(7ro )(i^o 
Moreover, using the estimate in Lemma 5.4, 

sup(e [4'Ul-M4'^))dfi 

nSN V J ^ ^ 



as n — )• oo. 



< sup 



nh. 



ntn+nhn 



n Jnt 



£; ( X^("'h^(„) ) ds ) -^ asc-^oo 



and 



E 



4')(l-V'c(4'^))d;Uo 



< E [XqI\Xo\>c\ ^ Q asc^oo. 



The last three displays imply the convergence in (8.14), and thus the result follows. 



D 
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Proof of Lemma 8.3. To show the tightness of {fJ-n}, it suffices to show that {I'n} is tight, 
where 

v^{A) := Ef,4A) = -— / P[{xil\,it - r}i")) G A]ds, A G B{S). 

"nrin Jnt„ 

However, the tightness of f„ is immediate in view of the tightness of {(X^". ,f/^".— r)s )}neN,seK+! 
which was proved in Proposition 5.1. 

Let ;U be a weak hmit point of /i„ and J : 5 — )■ M+ be defined by 

J{tt) := / e""[J(7r,ii) Aljdti, 
Jo 



J(7r,n):= sup ( |A(7r|^))| + |A^f | 



where 

J(7r,n) : = 

0<t<u 

Then J is continuous and bounded on S, and in order to show that fi is supported on Sq, it 
suffices to show that ^{J{-k) = 0) = 1 (see [5], p. 147). In turn, for the latter equahty to hold, 
it suffices to show that for all e > 0, Enn{J{n) > e) — )• 0, as n — )• oo. Now 

1 i-ntn+nhn /"OO 

i?/x„,(J(vr) >e) = — P{ e-"( sup [lAxf:;) | + \A{tJ, - r}W)|] Al)du> e)ds. 

^iT'n Jntn Jo 0<t<u 

(n) (n) 

Finally, noting that 77^^. — fis is continuous and using Lemma 5.6, we now have that the right 
hand side of the latter equation converges to as n — )• 00. The result follows. D 

Proof of Lemma 8.4. For a measure ly £ 1-^(3), let E'^ denote the expectation operator. For 
(a), we show that 

^m{^) ('/(vrj^))') _ ^m('^) ('/(vr^^))) = a.s. for ah bounded continuous / on S. (8.15) 

Note that 



^'^"(/(7riib)-^'^"(/(7r(i))) 



1 rnt„+nhn / ^ . ^ n \ 

— / E ^„J{X%.) - E ,.,f{xil\) ds 



nh. 

2t 

^ ^~ll/llsup -^ asn-^oo. 



This proves (8.15), and thus (a) follows. 



H An) 



Property (b) is immediate from the fact that fjl_^, — fjs is non-decreasing and continuous 
with initial value for each n. 

To prove (c), it suffices to show that for a.e. to and for every T,5 > 



E^^-^^£fsi4'^)d4'^^l)=0, 
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where fs is defined in (6.2). In turn, for the above equahty to hold, it suffices to show that 



T 



r(l)W-^(2) 



EiEn 1^ /5(7rW)dvr^^)Al 



0. 



The latter equality is immediate on noting that for every T,6 > 

i?(i?^"(^^/5(4i))d42)Af 

1 (-ntn+nhn / f-T 



nh 



■n Jnt 



E\ 1^ fs{xSu)d{v^r^.-vi^^){u)M]ds = 



and thus 



^M 



T 



fs{nP)dnPAl 



lim E 



T 



r(l)W-^(2) 



E^-i 1^ /,(7r^^^)dvr^^iAl 



Finally, consider (d). It suffices to show that for every < r < t < oo 



^'^(V'(^(i\vr(2))(0(vr«) - cPi4'^) - f £i0(7r(i))dn - <A'(l)[vrf - n^])) 



E 



where V' : "S" — >• M is an arbitrary bounded, continuous, Qs measurable map. 

Now fix such r, t, and '0. Assume without loss of generality that /x„ converges to /i. Combining 
this weak convergence with Lemma 5.5, we see that the left hand side of the last display is the 
limit of 



E 



£'^-(V^(7r(i),7r(2))(^(vrf )) - .^(^D) - / Ci4>{7r^^^)du - 0'(l)[7rf - 4^)])) 

Jr 

1 pnt„+nh„ 



nh. 



E 



n Jnt 



Ef,{^{xi':^.,f^: 



{«) An) _ '(n) 



m 



Hxi^t) - Hx^t) - CMxit)du-<P'im'si,-f,it] 



ds. 



To complete the proof, it suffices to show that the limit of the expression in the last display is 
0. Note that for </> e C^{[1, oo)) 

c^ixt^) - 4>ixir^) - f c^rU{xi-^)ds - v^r^mt^ 

Jo 



is a martingale, where P}"" (/>(!) := n[(f){l) — </)(l )] and 



'W, 



i(")^2„ 



£rV(x):=APn^xJ^ 

k=0 



X + 



k-l 



n 



(j){x) 



^,f\k). 



Thus, it suffices to prove that 

rntn+nhn 



lim 

n->-oo nh. 



1 r' 

hn Jntn 



E 



C^;'Uixit)-C^cPixit)]du 



ds = 
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and 



i-ntn+nhn 



lim 

n^oo nhn 



1 i-ntn+nrin , , 

T-\ pi"V(i)-0'(i: 



E 



.(n) _ .(n) 
"^s+t Vs+r 



ds = 0. 



The proofs for the last two equahties are completed as those for (6.3) and (6.4) upon using the 
uniform estimates in Corollary 5.1 and Lemma 5.5. D 



9. Appendix 



Proof of Proposition 2.1. We will only consider the case where {Xq^Yq) = {x,y) for some 
{x,y) G [l,oo) X [0, oo). The general case can be treated similarly. The unique solvability of 

(2.5) is an immediate consequence of the Lipschitz property of the Skorohod map, Lipschitz 
coefficients (note that f{x) = ^/x is a Lipschitz function on [l,oo)), and a standard Picard 
iteration scheme. 

We next argue the unique solvability of (2.6). For n G N, let a^"'' := inf{t > 0\Xt > n}, 
X^ := X^^^(„), and f^^'{y) := y V K Consider the equation 

y/") =Yo + C2A2 /* xi")/(")(yi"))d. + y^^ /* y#7w(^di?r . (9.1) 

JO Jo 

From the Lipschitz property of f^"^' and \//("^ it follows that, for each n, the above equation 
has a unique pathwise solution. Let r^*^) := inf{t > 0\Y^' = ^} and 9^^^ := r*^") A a'-'^\ Note 
that y^"-) solves (2.6) on [0,0'^"^]. Also, by unique solvability of (9.1), we have for all n G N, 
y("+i)(. A 6'(")) = y(")(- A 6'(")). Finally, letting 6l(°°) := lim„^oo ^l^"), the unique solution of 

(2.6) is given by the following: 



Ytiu;) 



'y}"^\uj), if < t < 6l(")(a;) for some n G N 

0, iu>e°^{uj). 



D 
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